Hyperbolic reflections as fundamental building blocks for multilayer optics 
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We reelaborate on the basic properties of lossless multilayers by using bilinear transformations. We 
study some interesting properties of the multilayer transfer function in the unit disk, showing that 
hyperbolic geometry turns out to be an essential tool for understanding multilayer action. We use a 
simple trace criterion to classify multilayers into three classes that represent rotations, translations, 
or parallel displacements. Moreover, we show that these three actions can be decomposed as a 
product of two reflections in hyperbolic lines. Therefore, we conclude that hyperbolic reflections 
can be considered as the basic pieces for a deeper understanding of multilayer optics. 



I. INTRODUCTION 

Although special relativity is perhaps the first theory 
that comes to mind when speaking about the interplay 
between physics and geometry, one cannot ignore that 
geometrical ideas are essential tools in the development 
of many branches of modern physics [l| . 

The optics of layered media is not an exception: in re- 
cent years many concepts of geometrical nature have been 
introduced to gain further insights into the behavior of 
multilayers. The algebraic basis for these developments 
is the fact that the transfer matrix associated with a loss- 
less multilayer is an element of the group SU(1, 1), which 
is locally isomorphic to the (2 + l)-dimensional Lorentz 
group S0(2, l).This leads to a natural and complete iden- 
tification between reflection and transmission coefficients 
and the parameters of the corresponding Lorentz trans- 
formation 

As soon as one realizes that SU(1, l)is also the ba- 
sic group of the hyperbolic geometry it is tempt- 
ing to look for an enriching geometrical interpretation 
of the multilayer optics. Accordingly, we have recently 
proposed 0, E| to view the action of any lossless mul- 
tilayer as a bilinear (or Mobius) transformation on the 
unit disk, obtained by stereographic projection of the 
unit hyperboloid of SO(2, 1). This kind of representa- 
tion has been previously discussed for, e.g., the Poincare 
sphere in polarization optics 0, 0, for Gaussian beam 
propagation in laser mode-locking and optical pulse 
transmission [lOj, and also in modelling visual process- 

ing mi- 

The point we wish to emphasize is that these bilin- 
ear transformations preserve hyperbolic distance between 
points on the unit disk. In Euclidean geometry any trans- 
formation of the plane that preserves distance can be 
written as a composition of reflections, which can be 
then considered as the most basic transformations. In 
fact, the composition of two reflections in straight lines 
is a rotation, or a translation, according these lines are 
intersecting or parallel. 

In hyperbolic geometry, each circle orthogonal to the 



boundary of the unit disk is a hyperbolic line and reflec- 
tions appear as inversions. However, we have an essential 
difference with the Euclidean case because there are three 
different kind of lines: intersecting, parallel, and ultra- 
parallel (which are neither intersecting nor parallel) Q. 
In consequence, the composition of two reflections in hy- 
perbolic lines is now a rotation, a parallel displacement, 
or a translation: these are precisely the transformations 
of the unit disk that preserve distance. 

A powerful way of characterizing transformations is 
through the study of the points that they leave invari- 
ant. For example, in Euclidean geometry a rotation can 
be characterized by having only one fixed point, while a 
translation has no invariant point. For a reflection the 
fixed points consist of all the points of a line (the reflec- 
tion axis). 

In this paper we shall consider the fixed points of 
the bilinear transformation induced by the multilayer, 
showing that they can be classified according to the 
trace of the multilayer matrix. From this viewpoint, 
the three transformations mentioned above; namely, ro- 
tations, parallel displacements, and translations appear 
linked to the fact that the trace of the multilayer transfer 
matrix has a magnitude lesser, equal or greater than 2. 

Since reflections appear as the basic building blocks 
of these geometric motions .12] 

, we show that any multi- 
layer action can be decomposed in terms of two inversions 
whose meaning is investigated. Such a decomposition is 
worked out for practical examples. This shows the power 
of the method and, at the same time, allows for a deeper 
understanding of layered media. 



II. MULTILAYERS AND THE UNIT DISK 

We first briefly summarize the essential ingredients of 
multilayer optics we shall need for our purposes 0. We 
deal with a stratified structure, illustrated in Fig. I, that 
consists of a stack of I, . . . , j, . . . , m, plane-parallel layers 
sandwiched between two semi-infinite ambient (a) and 
substrate (s) media, which we shall assume to be iden- 
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FIG. 1: Wave vectors of the input [E)^' and Ei '] and output 
[E^ and E { s +) ] fields in a multilayer sandwiched between two 
identical semi-infinite ambient and substrate media. 



tical, since this is the common experimental case. Here- 
after all the media are supposed to be lossless, linear, 
homogeneous, and isotropic. 

We consider an incident monochromatic linearly polar- 
ized plane wave from the ambient, which makes an angle 
#o with the normal to the first interface and has ampli- 
tude Ea . The electric field is either in the plane of 
incidence (p polarization) or perpendicular to the plane 
of incidence (s polarization). We consider as well another 
plane wave of the same frequency and polarization, and 

with amplitude El ' , incident from the substrate at the 
same angle [T^ . 

As a result of multiple reflections in all the interfaces, 
we have a backward-traveling plane wave in the ambi- 
ent, denoted Ea , and a forward- traveling plane wave 
in the substrate, denoted E^ + \ If we consider the field 
amplitudes as a vector of the form 



E 



E(- 



(1) 



which applies to both ambient and substrate media, then 
the amplitudes at each side of the multilayer are related 
by a 2 x 2 complex matrix M as , we shall call the multilayer 
transfer matrix, in the form 



The matrix M as can be shown to be Q 



Mas = 



1/T„ 

Ras /T as 



^asl as 





a 






0* 


a* 



(3) 



where the complex numbers R as and T as are, respec- 
tively, the overall reflection and transmission coeffi- 
cients for a wave incident from the ambient. Because 
I -Ras | 2 + I ^as | 2 = 1, we have the additional condition 
M 2 — \ft\ 2 = 1 ot j equivalently, det M as = +1 and then 
the set of lossless multilayer matrices reduces to the group 
SU(1, 1), whose elements depend on three independent 
real parameters. 

The identity matrix corresponds to T as = 1 and R as — 
0, so it represents an antireflection system. The matrix 
that describes the overall system obtained by putting two 
multilayers together is the product of the matrices repre- 
senting each one of them, taken in the appropriate order. 
So, two multilayers, which are inverse, when composed 
give an antireflection system |14| . 

In Refs. [2] and [3] we have proposed to view the multi- 
layer action in a relativisticlike framework. Without go- 
ing into details, it is convenient to characterize the state 
of the fields at each side of the multilayer by means of 
the "space-time" coordinates 



ha[E^*E^}, 
i[|£(+)| 2 -|£(->| 2 ], 



(4) 



for both ambient and substrate media. The coordinate e 3 
is the semi-difference of the fluxes (note that this number 
can take any real value) and, therefore, is constant be- 
cause the multilayer is lossless. In consequence, we have 
that 



( e °) 2 - (e 1 ) 2 - (e 2 ) 2 = (e 3 ) 2 = constant. 



(5) 



Equation (J5J defines a two-sheeted hyperboloid of radius 
e 3 , which without loss of generality will be taken hence- 
forth as unity 

A simple calculation shows that if one uses stereo- 
graphic projection taking the south pole (—1, 0, 0) as pro- 
jection center (see Fig. 2), the projection of the point 
(e , e , e 2 ) becomes in the complex plane 



E(- 



1 



E(+) 



(6) 



E„ 



E s 



(2) 



The upper sheet of the unit hyperboloid is projected into 
the unit disk, the lower sheet into the external region, 
while the infinity goes to the boundary of the unit disk. 

The geodesies in the hyperboloid are intersections with 
the hyperboloid of planes passing through the origin. 
Consequently, hyperbolic lines are obtained from these 
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defined as 




FIG. 2: Outline of the unit hyperboloid and a geodesic on 
it. We also show how a hyperbolic line is obtained in the 
unit disk via stereographic projection taking the south pole 
as projection center. 



by stereographic projection and they correspond to cir- 
cle arcs that orthogonally cut the boundary of the unit 
disk. 

It seems natural to consider the complex variables in 
Eq. © for both ambient and substrate. In consequence, 
Eq. J2Jl defines a transformation on the complex plane C, 
mapping the point z s into the point z a according to 



Za = $[M as ,Z s ] = 



/?* + a*z s 
a + f3z s ' 



which is a bilinear (or Mobius) transformation. The ac- 
tion of the multilayer can be seen as a function z a = f(z s ) 
that can be appropriately called the multilayer transfer 
function [(|. The action of the inverse matrix M^ 1 is 
z s = $[M~} , z a ). One can show that the unit disk, the 
external region and the boundary remain invariant under 
the multilayer action. 

For later purposes, we need the concept of distance in 
the unit disk. To this end, it is customary to define the 
cross ratio of four distinct points za, Zb, %c, and zd as 
the number 



{z A ,z B \z c ,z D ) = 



(z A - z c )/(z B - z c ) 
(z A - z D )/(z B - z D ) ' 



(8) 



which is real only when the four points lie on a circle or 
a straight line. In fact, bilinear transformations preserve 
this cross ratio [lfij . 

Let now z and z' be two points that arc joined by the 
hyperbolic line whose endpoints on the unit circle are 
E and E' . The hyperbolic distance between z and z' is 



d H (z,z') = -\hi(E,E'\z,z')\. 



(9) 



This can be seen as arising from the usual Minkowski 
distance in the unit hyperboloid (obtained through 
geodesies) by stereographic projection 0. The essential 
point for our purposes here is that bilinear transforma- 
tions are isometries; i.e., they preserve this distance. 



III. TRACE CRITERION FOR THE 
CLASSIFICATION OF MULTILAYERS 

Bilinear transformations constitute an important tool 
in many branches of physics. For example, in polarization 
optics they have been employed for a simple classification 
of polarizing devices by means of the concept of eigen- 
polarizations of the transfer function Q • The equivalent 
concept in multilayer optics can be stated as the field 
configurations such that z a — z s = zj in Eq. Q, that is 



z f = $[M aB ,z f ], 



whose solutions are 



i {-2z Im(a) ± v/[Tr(M as )p _ 4} . 



(10) 



(11) 



These values Zf are known as fixed points of the trans- 
formation <E>. The trace of M as provides then a suitable 
tool for the classification of multilayers 0] . 

When [Tr(M QS )] 2 < 4 the multilayer action is elliptic 
and it has only one fixed point inside the unit disk, while 
the other lies outside. Since in the Euclidean geometry 
a rotation is characterized for having only one invariant 
point, this multilayer action can be appropriately called 
a hyperbolic rotation. 

When [Tr(M as )] 2 > 4 the multilayer action is hyper- 
bolic and it has two fixed points both on the boundary 
of the unit disk. The hyperbolic line joining these two 
fixed points remains invariant and thus, by analogy with 
the Euclidean case, this action will be called a hyperbolic 
translation. 

Finally, when [Tr(M as )] 2 = 4 the multilayer action is 
parabolic and it has only one (double) fixed point on the 
boundary of the unit disk. This action will be called a 
parallel displacement. 

To proceed further let us note that by taking the con- 
jugate of M as with any matrix C <E SU(1, 1); i.e., 



= C M as C 



(12) 



we obtain another matrix of the same type, since 
Tr(M as ) = Tr(M as ). Conversely, if two multilayer ma- 
trices have the same trace, one can always find a matrix 
C satisfying Eq. l(T2"jl . 

The fixed points of M as are then the image by C of the 
fixed points of M QS . In consequence, given any multilayer 
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Rotation 
(elliptic) 



Translation 
(hyperbolic) 



Parallel displacement 
(parabolic) 




FIG. 3: Plot of typical orbits in the unit disk for: (a) canon- 
ical transfer matrices as given in Eq. 1131 and (b) arbitrary 
transfer matrices. 



matrix M as one can always reduce it to a 
of the following canonical forms: 



with 



K(<p) 
A(X) 



exp(i(p/2) 

exp(— iip/2) 

cosh(x/2) i sinh(x/2) 
—i sinh(x/2) cosh(x/2) 

l-iri/2 r)/2 
rj/2 l + irj/2 



(13) 



that have as fixed points the origin (elliptic), +i and 
— i (hyperbolic) and +i (parabolic), and whose physical 
significance has been studied before 0. The explicit 
construction of the family of matrices C is easy: it suffices 
to impose that C transforms the fixed points of M as into 
the ones of K(y>), A(x), or N(ry). 

The concept of orbit is especially appropriate for ob- 
taining an intuitive picture of these actions. We recall 
that given a point z, its orbit is the set of points z' ob- 
tained from z by the action of all the elements of the 
group. In Fig. 3. a we have plotted typical orbits for each 
one of the canonical forms K(ip), A(x), and For 
matrices K(ip) the orbits are circumferences centered at 
the origin and there are no invariant hyperbolic lines. For 
A(x), they are arcs of circumference going from the point 
+i to the point —i through z and they are known as hy- 
percicles. Every hypercicle is equidistant [in the sense of 
the distance |JjJJ|] from the imaginary axis, which remains 
invariant (in the Euclidean plane the locus of a point at 
a constant distance from a fixed line is a pair of paral- 
lel lines). Finally, for N(?7) the orbits are circumferences 
passing through the point +i and joining the points z 
and — z* and they are known as horocycles: they can be 
viewed as the locus of a point that is derived from the 
point +i by a continuous parallel displacement Q. 

For a general matrix M as the corresponding orbits can 
be obtained by transforming with the appropriate ma- 
trix C the orbits described before. In Fig. 3.b we have 
plotted typical examples of such orbits for elliptic, hy- 
perbolic, and parabolic actions. We stress that once the 
fixed points of the multilayer matrix are known, one can 



ensure that z a will lie in the orbit associated to z s . 

In the Euclidean plane any isometry is either a rota- 
tion, a translation, or a reflection. In any case, reflections 
are the ultimate building blocks, since any isometry can 
be expressed as the composition of reflections. In this Eu- 
clidean plane two distinct lines are either intersecting or 
parallel. Accordingly, the composition of two reflections 
in two intersecting lines forming an angle (p is a rotation 
of angle 2<p while the composition of two reflections in 
two parallel lines separated a distance d is a translation 
of value 2d. 

However, in the hyperbolic geometry induced in the 
unit disk, any two distinct lines are either intersecting 
(they cross in a point inside the unit disk), parallel (they 
meet at infinity; i.e., at a point on the boundary of the 
unit disk) , or ultraparallel (they have no common points) . 
A natural question arises: what is the composition of 
reflections in these three different kind of lines? To some 
extent, the answer could be expected: the composition 
is a rotation, a parallel displacement, or a translation, 
respectively. However, to gain further insights one needs 
to know how to deal with reflections in the unit disk. 
This is precisely the goal of next Section. 



IV. REFLECTIONS IN THE UNIT DISK 

In the Euclidean plane given any straight line and a 
point P which does not lie on the line, its reflected image 
P' is such that the line is equidistant from P and P' . In 
other words, a reflection is a special kind of isometry in 
which the invariant points consist of all the points on the 
line. 

The concept of hyperbolic reflection is completely anal- 
ogous: given the hyperbolic line V and a point P, to ob- 
tain its reflected image P' in T we must drop a hyperbolic 
line r from P perpendicular to T (such a hyperbolic line 
exists and it is unique) and extending an equal hyper- 
bolic distance [according to ©] on the opposite side of 
r from P. In the unit disk, this corresponds precisely to 
an inversion. 

To maintain this paper as self-containd as possible, let 
us first recall some facts about the concept of inversion. 
Let C be a circle with center w and radius R. An inver- 
sion on the circle C maps the point z into the point z' 
along the same radius in such a way that the product of 
distances from the center w satisfies 



\z' -w\ \z-w\ = R 2 , 
and hence one immediately gets 

R 2 R 2 + wz* - w*w 



z' = w + 



(14) 



(15) 



z* — w* z* — w* 

If the circle C is a hyperbolic line, it is orthogonal to the 
boundary of the unit disk and fulfills wio* = R 2 + 1 . In 
consequence 

, wz* — 1 



(16) 
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One can check 16] that inversion maps circles and lines 
into circles and lines, and transforms angles into equal 
angles (although reversing the orientation). If a circle C 
passes through the points P and P', inverse of P in the 
circle C, then C and C are perpendicular. Moreover, 
the hyperbolic distance JHJ) is invariant under inversions. 
This confirms that inversions are indeed reflections and 
so they appear as the most basic isometries of the unit 
disk. 

It will probe useful to introduce the conjugate bilinear 
transformation associated with a matrix M as as [compare 
with Eq. ©] 

Za = $*[M as ,z s ] = P*+ a ** : . (17) 
a + fjz* 

With this notation we can recast Eq. (|16|) as 

z' = $*[\ w ,z], (18) 

where the matrix \ w 6 SU(1, 1) associated to the inver- 
sion is 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



§ 2 (rad) 

FIG. 4: Plot of typical orbits in the unit disk for: (a) canon- 
ical transfer matrices as given in Eq. 1131 and (b) arbitrary 
transfer matrices. 



V. MULTILAYER ACTION AS COMPOSITION 
OF REFLECTIONS 



-i w*/R i/R 
—i/R i w/R 



(19) 



The composition law for inversions can be stated as fol- 
lows: if z' = z] and z" = then 



z" = $[\ w >\* w , z]. 



(20) 



To shed light on the physical meaning of the inver- 
sion, assume that incoming and outgoing fields are inter- 
changed in the basic configuration shown in Fig. 1. In 
our case, this is tantamount to reversing the time arrow. 
It is well know that given a forward-traveling field E^ + \ 
the conjugate field represents a backward phase- 

conjugate wave of the original field j2jj. In other words, 
the time-reversal operation can be viewed in this context 
as the transformation 



1 



(21) 



for both ambient and substrate variables; that is, it can 
be represented by an inversion in the unit circle. The 
transformed points lie outside the unit circle because, 
according to Eq. J3J this time reversal transforms the 
upper sheet into the lower sheet of the hyperboloid. 

Moreover, by direct inspection is easy to convince one- 
self that the matrix relating these time-reversed fields is 
precisely M* s and so the action can be put as 



(1/2-)* 



(3* +a*{l/z 8 )* 
a + p{l/z s )* ' 



(22) 



which expresses a general property of the time-reversal 
invariance in our model. 



As we have anticipated at the end of Section IIIII the 
composition of two reflections gives a rotation, a parallel 
displacement, or a translation, accordingly the two hy- 
perbolic lines are intersecting, parallel, or ultraparallcl, 
respectively. 

To illustrate how this geometrical scenario translates 
into practice, we consider an optical system versatile 
enough so as it could operate in the three regimes. To this 
end, we choose a Fabry-Perot-like system formed by two 
identical plates (each one of them with fixed phase thick- 
ness 5i) separated by a spacer of air with phase thickness 
62- This is a symmetric system for which R as and T as 
can be easily computed. By varying the spacer phase- 
thickness 82 we obtain the values of [Tr(M as )] 2 shown in 
Fig. 4. In all the examples we take as initial condition 
that in the substrate z s = 0.4exp(— in/ 3). 

First, we take 82 = 3 rad, so we are in the elliptic 
case ([Tr(M as )] 2 < 4). From the computed values of R as 
and T as one easily obtains the value z a = —0.4447 + 
0.4882i. The fixed point turns out to be z f = -0.3114 
and, in consequence, the multilayer action is a hyperbolic 
rotation around the center zt of angle 2<p, as indicated 
in Fig. 5. This multilayer action can be viewed as the 
composition of reflections in two hyperbolic lines Ti and 
T2 intersecting at Zf and forming an angle tp. The first 
inversion maps z s into the intermediate point Zj nt , which 
is mapped into z a by the second inversion. Note that 
there are infinity pairs of lines satisfying these conditions, 
but chosen arbitrarily one of them, the other is uniquely 
determined. Moreover, once these lines are known, they 
determine automatically the associated inversions. 

Next, we take the spacer phase thickness 82 — 1 rad, 
which corresponds to the hyperbolic case ([Tr(M as )] 2 > 
4), and we get z a = 0. 1567+0. 4464i. The fixed points are 
Zfx = -0.3695 - 0.9292i and z f2 = -0.3695 + 0.9292i. 



FIG. 5: Decomposition of the multilayer action in terms of 
two reflections in two intersecting lines for the same multilayer 
as in Fig. 4 with 82 = 3 rad (elliptic case). 




FIG. 6: Decomposition of the multilayer action in terms of 
two reflections in two ultraparallel lines for the same multi- 
layer as in Fig. 4 with 82 = 1 rad (hyperbolic case). 

There are no invariant points in the unit disk, but the 
hyperbolic line joining z/i and z/ 2 is the axis of the 
hyperbolic translation. In Fig. 6, we have plotted the 
hypercicle passing through z a and z s . The multilayer 
action can be now interpreted as the composition of two 
reflections in two ultraparallel hyperbolic lines Ti and T 2 
orthogonal to the translation axis. If Ti and r 2 intersect 
the hypercicle at the points z\ and z 2 , they must fulfill 

dn{z a , z s ) = 2da(zi, z 2 ), (23) 

in complete analogy with what happens in the Euclidean 
plane. Once again, there are infinity pairs of lines fulfill- 
ing this condition. 



FIG. 7: Decomposition of the multilayer action in terms of 
two reflections in two parallel lines for the same multilayer as 
in Fig. 4 with 82 = 0.4328 rad (parabolic case). 



Finally, we take 82 — 0.4328 rad, so we are in the 
parabolic case ([Tr(M as )] 2 = 4), and z a = -0.1615 + 
0.6220i. The (double) fixed point is z f = -1. In Fig.7 
we have plotted the horocyle connecting z s and z a and 
the fixed point. Now, we have the composition of two 
reflections in two parallel lines Ti and r 2 that intersect 
at the fixed point Zf and with the same constraint ill'MI) 
as before. 



VI. CONCLUDING REMARKS 

In this paper, we have proved a geometric scenario to 
deal with multilayer optics. More specifically, we have re- 
duced the action any lossless multilayer (no matter how 
complicated it might be) to a rotation, a parallel dis- 
placement or a translation, according to the magnitude 
of its trace. These are the basic isometries of the unit disk 
and we have expressed them as the composition of two 
reflections in intersecting, ultraparallel, or parallel lines. 
There is no subsequent factorization in simpler terms so, 
reflections are the most basic motions one can find in the 
unit disk. 

We hope that this approach will complement the more 
standard algebraic method in terms of transfer matrices, 
and together they will aid to obtain a better physical and 
geometrical feeling for the properties of multilayers. 

Finally, we stress that the benefit of this formulation 
lies not in any inherent advantage in terms of efficiency in 
solving problems in layered structures. Rather, we expect 
that the formalism presented here could provide a general 
and unifying tool to analyze multilayer performance in a 
way closely related to other fields of physics, which seems 
to be more than a curiosity. 
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